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The dynamics of fronts, such as chemical reaction fronts, propagating in two-dimensional fluid 
flows can be remarkably rich and varied. For time-invariant flows, the front dynamics may simplify, 
settling in to a steady state in which the reacted domain is static, and the front appears “frozen”. 
Our central result is that these frozen fronts in the two-dimensional fluid are composed of seg¬ 
ments of burning invariant manifolds —invariant manifolds of front-element dynamics in xyO- space, 
where 0 is the front orientation. Burning invariant manifolds (BIMs) have been identified previously 
as important local barriers to front propagation in fluid flows. The relevance of BIMs for frozen 
fronts rests in their ability, under appropriate conditions, to form global barriers, separating reacted 
domains from nonreacted domains for all time. The second main result of this paper is an under¬ 
standing of bifurcations that lead from a nonfrozen state to a frozen state, as well as bifurcations 
that change the topological structure of the frozen front. Though the primary results of this study 
apply to general fluid flows, our analysis focuses on a chain of vortices in a channel flow with an 
imposed wind. For this system, we present both experimental and numerical studies that support 
the theoretical analysis developed here. 

PACS numbers: 47.70.Fw, 82.40.Ck, 47.10.Fg 


I. INTRODUCTION: REACTING FLOWS AND 
FROZEN FRONTS 

The evolution of an autocatalytic reaction A-\-B -A 2 A 
in a spatially-extended system is characterized by the 
propagation of reaction fronts that separate the species A 
and B. The motion of these fronts is well-understood for 
react ion-diffusion systems in the absence of any substrate 
flow. The effects of fluid motion on fronts in the more 
general advection-reaction-diffusion system have only re¬ 
cently received significant attention. This is somewhat 
surprising, given the applicability of advection-reaction- 
diffusion to a wide range of systems, including microflu¬ 
idic chemical reactors [1], plasmas [2], the dynamics of 
ecosystems in the oceans (e.g., plankton blooms) [3], 
cellular- and embryonic-scale biological processes [4, 5], 
and the propagation of diseases in society [6] . It has been 
recently proposed that the motion of reaction fronts in 
fluid flows may be dominated by the presence of burning 
invariant manifolds (BIMs), which act as one-way barri¬ 
ers to advancing fronts [7, 8] . The existence of BIMs and 
their function as one-way barriers has been verified exper¬ 
imentally in time-independent and time-periodic vortex 
chain flows, as well as 2D disordered vortex flows [9]. 

Experiments have shown that reaction fronts tend 
to pin to vortex structures in the presence of an im¬ 
posed wind [10]. These fronts neither propagate forward 
against the wind nor are blown backwards, but remain 
“frozen”. This behavior is surprisingly robust, occurring 
over more than an order of magnitude of wind speeds and 
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FIG. 1. Sequence showing the evolution of a triggered re¬ 
action front. The maximum fluid vortex speed (in the ab¬ 
sence of wind) is U — 1.4 mm/s, and the wind speed is 
V w = 0.90 mm/s. The images in the sequence are separated 
by 10 s. 


a variety of underlying flows ranging from confined vor¬ 
tex chains to extended, spatially-random flows. Figure 1 
shows a sequence of images from experiments showing the 
evolution of a triggered, autocatalytic reaction front in a 
vortex chain with wind. The front eventually stabilizes 
and remains fixed for the duration of the experiment. 

In this paper, we use the theory of BIMs to characterize 
these frozen fronts (FFs). FFs occur when a BIM spans 
the entire width of the system with no changes in block¬ 
ing direction, or when there is a combination of over¬ 
lapping BIMs, with the same blocking directions, that 
together span the system. In either of these situations, 
the shape of the FF is determined by the shape of the 
BIMs responsible. We illustrate the creation of FFs and 
changes in their structure by increasing the wind applied 
to a canonical base flow (the alternating vortex chain) 
with a propagating chemical reaction. We present both 
experimental and numerical studies of this system. 

This paper is organized as follows. We begin in Sec. II 
by presenting experiments involving reaction fronts in a 
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particular quasi-two-dimensional fluid flow—the “windy 
alternating vortex chain”. The images in this section il¬ 
lustrate the behavior of FFs under an imposed wind of 
various strengths. Section III recalls some basic aspects 
of burning invariant manifolds (BIMs)—geometric struc¬ 
tures that govern the progress of fronts in fluid flows— 
including the three-dimensional dynamics of front ele¬ 
ments and fixed points of this system. Next, Sec. IV 
connects the previous two sections by showing that FFs 
are composed of BIM segments. Section V considers FFs 
in a numerical model of the experimental flow. It paral¬ 
lels Sec. II by increasing the applied wind and observing 
the resulting changes in the FFs. Here we discuss in 
detail the various FF topologies and the dynamical sys¬ 
tems mechanism underlying the transitions which con¬ 
nect them. There are four appendices. Appendix A in¬ 
troduces a two-dimensional invariant surface of “sliding 
fronts”, which is used to prove several key results in the 
paper. Appendix B establishes the stability condition 
that frozen fronts must satisfy. Appendix C provides 
some technical analysis concerning the structure of frozen 
fronts at infinity. Finally, Appendix D examines attract¬ 
ing fixed points for this system. 


II. EXPERIMENTS: WINDY ALTERNATING 
VORTEX CHAIN FLOW 


The alternating vortex chain fluid flow has been the 
subject of much study, both theoretical and experimen¬ 
tal. It has been used as a model of a two-dimensional 
cross-section of Rayleigh-Benard (thermal) convection 
[11-14] and Taylor-Couette vortices [15], and can be used 
to model vortex chains and streets in oceanic and at¬ 
mospheric flows [16, 17]. The alternating vortex chain 
has been used to study enhancement of long-range, fluid 
transport in cellular flows [18-21]. More recently, it has 
been used repeatedly in studies of chemical front propa¬ 
gation in advection-reaction-diffusion systems [7, 22-26]. 
Here we modify this flow by adding a uniform “wind”, 
creating the windy alternating vortex chain [10]. 



FIG. 2. Experimental apparatus, (a) Exploded view of alter¬ 
nating fluid vortices above array of magnets. Current through 
fluid induces Lorenz force, (b) Side view of the apparatus. 
A 2 mm thick layer of an electrolytic fluid is contained in 
an acrylic box. The chain of Nd-Fe-Bo magnets moves on a 
translation stage below the box. 


A. Experimental setup 

The experiments are conducted in a quasi-2D flow 
composed of a chain of vortices in a thin (2 mm) fluid 
layer. The flow is produced using a magnetohydro dy¬ 
namic forcing technique, as shown in Fig. 2. A chain 
of permanent 1.9cm-diameter Nd-Fe-Bo magnets sits be¬ 
low the fluid layer, thereby imposing a spatially-varying 
magnetic field. An electric current is passed though this 
electrolytic fluid, generating Lorentz forces on the fluid. 
In conjunction with rigid, plastic side-walls that bound 
the region of interest, the result is an alternating chain of 
well-controlled vortices. The magnets are mounted on a 
translation stage; motion of the translation stage results 
in motion of the magnets and, consequently, the fluid vor¬ 
tices. In these experiments, we move the magnets (and 
the vortices) with a constant speed V w . In the reference 
frame moving with the magnets, the flow is a stationary 
chain of vortices with an imposed, uniform wind of speed 

V w . 

The fronts are produced in the experiments with the 
excitable, ferroin-catalyzed Belousov-Zhabotinsky chem¬ 
ical reaction [27, 28]. At the beginning of an experimental 
run, the ferroin indicator in the solution is in its reduced 
(orange) state. A reaction is then triggered by briefly 
dipping a silver wire into the fluid. The silver oxidizes 
the ferroin in its vicinity, changing the local indicator to a 
blue-green color. The oxidized indicator in turn oxidizes 
the ferroin of its neighbors, resulting in a blue-green reac¬ 
tion front that steadily propagates outward from the trig¬ 
ger point with a roughly constant propagation speed Vo. 
For all experiments presented in this article, the propa¬ 
gation speed is Vo = 0.07 mm/s. The front is a pulse-like 
front—behind the leading edge of the front, the reaction 
relaxes back to its reduced (orange) state and can be re¬ 
triggered. Previous studies [24-26] have indicated that 
the behavior of the leading edge of these pulse-like fronts 
in a fluid flow is identical to the behavior of the leading 
edge of a burn-type reaction. (Burn-type reactions do 
not relax back, rather A + B — 2A and stays that way.) 

B. Experimental results 

We focus on the behavior of the leading edge of the re¬ 
action front that propagates against the imposed wind. 
(In the lab frame, these fronts propagate in the direction 
of the imposed motion of the vortex cores.) An example 
of a typical experiment is shown in Fig. 3. As viewed in 
the laboratory reference frame (Fig. 3a), the front con¬ 
tinually propagates in both directions; in the reference 
frame moving with the vortices (Fig. 3b) the right-most 
edge of the reaction front converges to a steady-state sta¬ 
tionary shape that remains fixed for the duration of the 
run. From here on, we use the expression “wind” V w to 
refer to either the translational speed of the vortices in 
the lab frame or the speed of the uniform wind in the 
vortex reference frame. 
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FIG. 3. Sequences showing the evolution of a reaction front 
in a vortex chain, (a) Lab frame, with the vortices moving to 
the right, (b) Reference frame moving with the vortices. In 
this frame, the vortices are stationary and a wind blows across 
the vortices toward the left. U = 1.4 mm/s, V w — 0.30 mm/s. 
Images in the sequences are separated by 20 s. 


The propagation of a reaction front in the alternating 
vortex flow in the absence of an imposed wind has been 
discussed in detail in previous papers [25, 26, 29, 30]. 
The reaction front is carried around each vortex with the 
flow and “burns” across the separatrix from one vortex 
to the next, resulting in long-range propagation that is 
significantly faster than the reaction-diffusion speed Vo 
in a static fluid. The long-term average front speed is 
independent of the initial stimulation. 

If a uniform wind V w < Vo is applied (i.e., the wind 
speed is smaller than the reaction-diffusion speed), the 
reaction front still propagates to the right against the 
wind, although the long-range propagation speed is re¬ 
duced. At V w = Vo, there is a transition where the front 
neither advances against the wind nor is blown backwards 
[10]. Figure 4 shows a sequence for a reaction front trig¬ 
gered in a flow with wind V w just below Vo (Fig. 4a-4c) 
and V w just above Vo (Fig. 4d-4f). The shape of the FF 
is not arbitrary; rather a wide range of initial stimula¬ 
tions will result in fronts that converge onto the same 
structure. For V w = Vo, the shape of the FF corresponds 
well with the advective separatrix having V w = 0. 

As the strength of the imposed wind is increased, the 
shape of the FF evolves considerably. Figure 5 shows 
time-averaged images of the steady-state reaction fronts 
for several different wind speeds. With increasing wind 
speed, the contact point of the FF with the upper bound¬ 
ary does not move much. There is also a shift-flip sym¬ 
metry apparent in Fig. 5b-5g; for every FF originating 
from a contact point there is a flipped version of the 
same structure originating from a contact point one vor¬ 
tex width leftward. Consequently, for any wind speed, 
the leading edge of the front could be pinned to any one 
of these contact points; i.e., any FF could be replaced by 
the same shape, shifted by one vortex and flipped verti¬ 
cally. 

The shift-flip symmetry is also relevant to a change in 
the structure of the FF as the wind speed is increased. 



FIG. 4. Two sequences demonstrate front evolution near crit¬ 
ical wind speed. The maximum fluid vortex speed (in the ab¬ 
sence of wind) is U — 0.7 mm/s. For wind value V w = 0.085 
mm/s, the front is (a) nearly vertical, (b) finds a small passage 
into the next right vortex, and (c) fills in the right vortex con¬ 
tinuing down the channel. For wind value V w = 0.090 mm/s, 
the front is (d) nearly vertical, (e) does not find passage to 
the right, and (f) remains unchanged from (e)—a frozen front. 
The time between images is 40 s in both sequences. Note that 
there is a small amount of experimental noise that increases 
the transition slightly above V w = Vo = 0.07 mm/s. 


The front develops a point, or corner, with an apparently 
discontinuous derivative (Fig. 5b). This point moves left¬ 
ward for larger and larger wind speeds (Fig. 5c-5g). This 
concave corner first appears near the downwind contact 
point (one vortex width downwind in Fig. 5b) and moves 
away from the channel wall. In this situation, the FF is 
composed of a combination of smooth curves that origi¬ 
nate at different contact points. 

Above a minimum wind speed, the shape of the FF 
is no longer uniquely determined (modulo the flip-shift 
symmetry); rather, more than one front shape is possible, 
depending on the manner in which the front is triggered 
(Fig. 6). It is possible to trigger a reaction front that pins 
only to the structure emanating from a single contact 
point, as in Fig. 6a. But the same flow allows for other 
FFs, such as in Fig. 6b. The number of different possible 
FF shapes increases with the wind. As can be seen in 
both Fig. 5 and Fig. 6, the front shapes are stretched 
out significantly with increasing wind speed, spanning 
more and more vortex cells. For all except the smallest 
wind speeds, a FF can be composed of structures pinned 
onto adjacent vortex contact points, as in Figs. 5b-5g 
and Figs. 6b and 6c. For larger wind speeds, additional 
FF shapes are possible. As an example, Fig. 6d shows a 
FF composed of two structures originating from contact 
points separated by 5 vortex widths. 

Experimentally, the more complex steady-state front 
shapes are often found by simultaneously triggering the 
reaction in multiple locations. However, these complex 
shapes appear to be sometimes accessible with even a 
single, well-placed trigger. A more detailed theoretical 
treatment of these “basins of attraction” is in prepara¬ 
tion. 

For large enough wind, the stable state is lost com¬ 
pletely, with the front being “blown backwards” down- 
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FIG. 5. (Color online.) Time-averaged images of steady- 
state reactions for several wind speeds. U = 1.4 mm/s for 
all. V w = (a) 0.15 mm/s, (b) 0.16 mm/s, (c) 0.20 mm/s, (d) 
0.30 mm/s, (e) 0.60 mm/s, (f) 0.90 mm/s, and (g) 1.2 mm/s. 
Arrows indicate apparent discontinuities in the FF tangent 
direction. 



FIG. 6. Multiple FFs are realized with the same flow, de¬ 
pending on how the reaction is triggered. U = 1.4 mm/s for 
all. V w is 0.90 mm/s for (a,b) and 1.2 mm/s for (c,d). 


wind. A complete parameter space showing the range of 
wind speeds for FFs can be found in Ref. [10]. 


III. BIM REVIEW 

We model advection-reaction-diffusion systems, such 
as the above experiments, by considering only the front. 
This amounts to taking the so-called “sharp-front”, or 
geometric-optics limit. While some other studies have 
made use of a grid-based computational scheme [23, 24], 


focusing on the front is numerically economical and the¬ 
oretically insightful. By assuming that the front pro¬ 
gresses in a curvature-independent way [31], the front 
may be regarded as the collection of independent front 
elements that comprise it. Although not crucial to 
the basic ideas here, we also assume that the “burning 
speed” [32] (i.e. front propagation speed in the local fluid 
frame) is homogenous and isotropic. 

A front is the oriented boundary of a burned region 
with orientation defined by the normal vector n pointing 
away from the burned region. (We can also refer to the 
orientation using the tangent vector g where n x g = +1, 
i.e. pointing out of the plane.) Denoting by r the ap¬ 
position of a front element and by 6 the angle from the 
x-axis to g, a front is a curve in xyO -space that satisfies 
the front-compatibility criterion , 

!“«*>■ (i > 

where A is some smooth parameterization of the curve. 
The above assumptions lead to the following three- 
dimensional ODE governing the evolution of an individ¬ 
ual front element (r (t),6(t)). 


r = u + r 0 n, (2a) 

6 = - hiUijgj , (2b) 

where u is the prescribed fluid velocity field, which is 
nondimensionalized by dividing by t/, the maximum fluid 
vortex speed in the absence of wind. That is, in the ab¬ 
sence of wind, the maximum value of u is unity. Here, 
vq = Vo/U is the nondimensionalized front propagation 
speed in the comoving fluid frame. The position variable 
r is scaled so that the width of each vortex and of the 
channel is unity. Time is scaled by the advection time 
D/C/, where D is the (dimensionful) vortex width. Note 
that g = (cos #, sin 0) and n = (sin 6, — cos 6) indicate 
the tangent to the front element and the normal direc¬ 
tion (propagation direction), respectively. Furthermore, 
Uij = dui/drj and repeated indices are summed. The 
total translational motion of a front element is the vector 
sum of the fluid velocity and the front propagation veloc¬ 
ity in the fluid frame, Eq. (2a). The change in orientation 
is determined entirely kinematically; Eq. (2b) describes 
the angular velocity of a material line embedded in the 
fluid. 

Invariant manifolds of the full 3D ( xyO ) dynamics, 
Eq. (2), depend upon both the fluid flow and front prop¬ 
agation, and therefore differ from the invariant mani¬ 
folds of the underlying advection dynamics. We focus 
on the ID unstable manifolds attached to the burning 
fixed points (BFPs)—i.e. fixed points of Eq. (2)— that 
are of stability type stable-stable-unstable (SSU). We 
call these burning invariant manifolds (BIMs). It has 
been demonstrated theoretically and experimentally that 
these BIMs are “one-way” barriers to front propagation 
in flows (Fig. 7). That is, they prevent reactions from 
crossing in one direction but allow them to cross in the 
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FIG. 7. (Color online.) Evolution of reaction front (blue 
to green) in two counter-rotating vortices. Stimulation on 
lower left grows while being acted on by the flow. Two BIMs 
(red) emanate from BFPs on the bottom channel wall. The 
“burning direction” of each BIM is indicated by red triangles. 
The reaction passes through oppositely oriented BIM, but is 
blocked by cooriented BIM. Finally the reaction front wraps 
around cusp of right BIM. 


other. It is somewhat surprising that these codimension- 
two manifolds are in fact barriers. BIMs are not generic 
curves through xyO- space; they obey the front compati¬ 
bility criterion 1 [8]. All fronts, including BIMs, obey the 
front no-passing lemma: no front can overtake another 
front from behind. 

An interesting consequence of the front propagation 
dynamics is the ability to create cusps in fronts and in the 
BIMs. In time-independent flows, cusps mark a change 
in the bounding nature of BIMs. Figure. 7 illustrates 
the evolution of a small circular front (lower left, blue). 
During its evolution (blue to green), it passes through the 
left BIM (red) because of their opposite orientation. It 
then presses up against the right BIM (cooriented) and 
follows closely until reaching the BIM cusp where the 
BIM’s relative orientation changes, thus allowing passage 
of the reaction front. We define the BIM core as the 
BIM segment that includes the BFP and extends in both 
directions until reaching either a cusp, a new BFP, or 
infinity. 


IV. FROZEN FRONTS: BASIC THEORY 

Consider a fluid domain D that is connected, but not 
necessarily simply connected. In this paper, we focus on 
a channel flow where D = M ® [0,1], but the results ob¬ 
tained in this section are general. We now introduce a 
more precise mathematical definition of frozen front than 
the more intuitive definition used thus far. First, we de¬ 
fine frozen domain as a burned subdomain of D that 
is invariant under the burning dynamics and stable to 
perturbation [33]. (See App. B for a precise discussion of 
this notion of stability). Since the fluid is incompressible, 
neither the frozen domain nor its complement may be of 
finite area. A frozen front (FF) is the oriented bound¬ 
ary of a frozen domain that separates the burned from 




allowed transitions 


FIG. 8. An arbitrary front element exists in one of four states 
with respect to a burned region: A, unburned region; B, on 
the boundary of the burned region with non-outward-normal 
burning direction; C, on the boundary and oriented in the 
outward normal direction; D, inside the burned region. The 
diagram on the right indicates how the state of a front element 
may change as it coevolves with the burned region. These 
same dynamics hold between a front element and the fluid 
domain boundary. 


the unburned fluid. (The frozen domain boundary that 
coincides with the boundary of D, i.e. a domain wall, is 
then not considered part of the frozen front.) As with 
any front, we choose the orientation of the FF to be a 
unit vector normal to the FF pointing outward from the 
burned region. Since the frozen domain is unbounded, 
the FF cannot be a closed curve. 

Consider a particular FF F as a curve in xyO- space. An 
individual front element on F can evolve into the interior 
of the frozen domain, but not vice versa (Fig. 8). Since 
the frozen domain is invariant, the time evolution of F 
under Eq. (2) includes F for any time t. In other words, 
the backward trajectory of any point on F remains on 
F. Thus the FF must be the union of segments of front 
element trajectories, and is hence a piece-wise smooth 
curve. Each segment follows a trajectory from t = — oc 
to some t = tf. This implies each segment lies within the 
unstable manifold emanating from a fixed point, which 
may be at infinity (see App. C). 

On a smooth segment of FF a front element is either 
a fixed point of the flow, or it “slides” along the seg¬ 
ment satisfying rocg. Any FF can thus be decomposed 
into a collection of these sliding fronts (App. A). Here 
we summarize the geometry of sliding fronts detailed in 
App. A. First, sliding fronts only exist in the domain 
where |u| > v$. We refer to this domain as the fast zone 
FZ, and the complementary domain as the slow zone SZ. 
In the FZ, the structure of the sliding fronts can be sim¬ 
ply characterized. At every point in the FZ interior, there 
are two allowed sliding front orientations characterized by 
the angle, 


P = arccos(—no/|u|), (3) 

between the front propagation direction n and the fluid 
flow u (Fig. 9, Lemma 8). In the limit uo/|u| —^ 0, the 
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FIG. 9. A generic intersection of two sliding fronts (one black, 
one gray). Each sliding front’s propagation vector aoh cancels 
the normal component of the fluid velocity u, leaving only 
motion tangent to the front. The two orientations (black and 
gray) are symmetric about u. 


two sliding fronts become parallel (burning in opposite 
directions) and align with the streamlines, thus recover¬ 
ing the advective case. We refer to these two choices of 
orientation as and ” corresponding to sgn(r • g). 
Each choice of orientation defines a set of sliding fronts 
whose projection foliates the FZ. When the sliding fronts 
are considered as curves in xyO- space, they foliate a two- 
dimensional surface which is a double-branched covering 
of the FZ. (See Figs. 23 and 24 for examples.) 

Consider a burned region bounded by two sliding fronts 
(on different branches) that meet at a point as in Fig. 9. 
In principle, the burned region may be either locally con¬ 
cave or locally convex at this point. However, the convex 
case is not relevant to FFs because any convex corner 
will be smoothed out after an arbitrarily short evolu¬ 
tion. Therefore, in the interior of the FZ, a FF is simply 
a union of smooth curves that meet at concave angles 
specified by the local burning-to-fluid-speed ratio vq/\u\. 
In the limit |u| vq, the two branches meet on the 
boundary of the SZ. At all such points of the boundary, 
two sliding fronts meet with burning directions n aligned. 
There are two cases to consider. 

In the first case, assume n is not perpendicular to the 
SZ. Then the sliding front trajectory passes through the 
fold joining the two branches in such a way that it forms 
a cusp in the xy -plane (Fig. 10). We observed above that 
cusps mark a change in the bounding behavior of BIMs. 
This change occurs at cusps along any sliding front (in¬ 
cluding BIMs), which implies that a FF cannot contain 
a cusp. Figure 11 illustrates why; it shows the two pos¬ 
sible burned regions that would be bounded by such a 
cusp. In both cases, one segment of the sliding front has 
a burning direction incompatible with, i.e. pointing into, 
the proposed burned region. 

Referring to Fig. 12, as n becomes perpendicular to the 
SZ at the point x, the cusp becomes tangent to the SZ. 
By symmetry, a cusp also approaches x from the other 
side. 

In the second case, where n is perpendicular to the SZ, 
Ref. [8] showed that the sliding front must meet the SZ 
at a BFP x. This could be thought of as the meeting of 
two cusps (Fig. 12). Each segment of the cusp on the left 



FIG. 10. (Color online.) The BIM (red) is a smooth curve 
in xyO- space. Its projection (green) onto the xy-pl&ne has a 
cusp on the boundary of the SZ. 



FIG. 11. (Color online.) A sliding front (red) with a cusp 
cannot bound a burned region (gray). Either choice of shading 
leads to an incompatibility in front orientation in which one 
piece of the front points into the burned region. 


pairs with its symmetric segment on the right to form 
a smooth curve in xy -space passing through x. Each 
of these two combined curves is a ID stable or unstable 
manifold of x. There are four possible stability types of 
BFPs in xyO- space: SSS, SSU, SUU, and UUU. These 
are illustrated in Fig. 25. For SSU and SUU BFPs, the 
dynamics restricted to the sliding surface is of stability 
SU (Lemma 1). Figure 25 illustrates the ID stable and 
unstable manifolds attached to such BFPs. For SSS and 
UUU points, the dynamics within the constraint surface 
is of stability SS and UU respectively (Lemma 1). Since 
the BFP is either a sink or source in this case, it is met 
by an infinite number of sliding trajectories. 

Only two of the four stability types can occur on a 
frozen front. Suppose a frozen front is tangent to a SZ 
at a BFP where the burning direction is into the SZ, as 
for SUU or UUU stability types. Though the burned re¬ 
gion behind the BFP does not intersect the SZ, a small 
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FIG. 12. (Color online.) Sliding fronts (black) strike the 
SZ, forming cusps on either side of an SSU BFP x. As the 
cusps on either side approach x, they become more horizontal, 
eventually joining tangent to each other at x. 

perturbation of the burned region at the BFP can inter¬ 
sect the SZ. Once any of the SZ is burned, the entirety 
of the SZ must eventually be burned and remain burned 
forever (Lemma 2). Since we require frozen fronts to 
be stable under small perturbations (App. B), SUU and 
UUU BFPs cannot occur on a FF. 

The two remaining stability types SSU and SSS can 
exist on a FF. We previously showed that the FF con¬ 
sists of unstable manifolds. Only the SSU points have 
unstable manifolds. Finally, since cusps are not allowed 
on FFs (shown earlier), we have one of the main results 
of this paper. 

Proposition 1. Frozen fronts are built from BIM cores. 
More precisely, each frozen front is generated by some 
set Sff of SSU BFPs. The frozen front is obtained by 
tracing the unstable manifold from each point in Sff 
until one of three things occurs: it intersects any other 
BIM core emanating from Sff / it intersects any domain 
boundary; or it terminates at an SSS BFP. 

So far we have focused our attention on the generation 
of the FF from BFPs. Here we shift our attention to 
consider how the sliding segments of a FF end. We have 
already discussed the most common case where segments 
intersect at a concave corner. There exist two other pos¬ 
sibilities, termination on an SSU or SSS BFP. 

An SSU BFP has a stable manifold that contains an 
incoming sliding front. A FF can therefore contain a 
segment which is a heteroclinic connection consisting of 
a sliding front between SSU points. Figure 13b shows 
two SSU BFPs joined by such a connection flowing from 
the upper to the lower BFP. This configuration is a FF; 
in particular, it is stable to perturbations of the burned 
region (App. B). In one sense, the FF is also structurally 
stable because generic perturbations of the flow yield 
frozen domains with a similar shape (Figs. 13a and 13c). 
In another sense, it is not structurally stable, because 
generic perturbations break the heteroclinic connection, 
thus altering the dynamics along the front. Some of these 
perturbations cause the lower SSU BFP to fall behind 
the FF (Fig. 13a), while other perturbations cause it to 



FIG. 13. (Color online.) The SSU-SSU connection is not 
structurally stable as the wind speed is varied. Nevertheless, 
the frozen domain (light gray) varies continuously. SZs are 
dark gray, (a) The relation v w < v c places the lower SSU 
point behind the FF attached to the upper SSU point, (b) 
The equality v w — v c makes the unstable manifold from the 
upper SSU point coincide with the stable manifold of the lower 
point, (c) The relation v w > v c pushes the lower SSU point 
ahead, placing it and its BIM on the FF. The FF is now 
composed of two BIMs meeting at a concave corner. 

push through, and in doing so contribute a segment of 
unstable sliding front to the FF (Fig. 13c). As seen in 
Figs. 13b and 13c, both of these perturbations return the 
system to the generic case. So while SSU BFPs can exist 
as “termination points” along a FF, this is not generic. 

Finally, we consider the SSU to SSS connection. The 
SSS point attracts all points within a 3D neighborhood 
and, therefore, it attracts all sliding fronts within some 
neighborhood on the invariant sliding surface. It might 
then seem that this SSS point can be on a FF containing 
any of these incoming sliding fronts. However, the sliding 
front must reach the SSS point without having formed a 
cusp. This can only happen if the eigenvalues of the SSS 
point are real (see App. D). Such SSS points do exist, 
albeit for what appears to be a small parameter range. 

V. THEORY: WINDY ALTERNATING VORTEX 
CHAIN FLOW 

We continue our discussion of FFs using a simple nu¬ 
merical model of the experimental fluid flow. 

A. Numerical model 

The stream function that describes the flow is 

T = — sin(7nr) sin(7 xy) — v w y , (4) 

7T 

where u x = d^/dy and u y = —d^t/dx. This model 
has been used in several previous studies, on both fluid 
mixing and reacting flows, yielding reasonable agreement 
with experiment. Our intent here is to illustrate the the¬ 
ory of frozen fronts for a particular fluid flow, and to re¬ 
produce basic features of the experimental flow in Sec. II. 

There is a weak three-dimensional component to the 
vortex flow due to Ekman pumping that carries fluid to- 

















ward the vortex centers at the bottom of the fluid layer 
and up through the vortex cores [34]. This effect is not 
included in the model. Also, while the model has free- 
slip boundary conditions, this is certainly not true in the 
experiment. Nevertheless, the simplified free-slip model 
of Eq. (4) has been used successfully in modelling several 
experiments on passive transport and front propagation 
in vortex flows [14, 20, 22, 24, 25]. The last term in 
Eq. (4) numerically models the fluid wind observed in 
the moving frame of the vortices. 


B. Dynamical systems analysis 

We begin by considering a flow where the wind is of 
insufficient strength to produce a FF (Fig. 14a). The 
streamlines indicate that this is essentially a vortex flow, 
but with a sinuous, left-moving jet. In Fig. 14b a small 
circular stimulation (purple) is made in the lower left. 
This circle evolves outward to the left and right while 
being deformed by the flow. The rightward moving front 
is able to make slow progress “upwind”. Notice though 
that it is blocked at the vortex boundary near the bottom 
and middle, and must wind around through the top of 
the channel. In this way, the reaction continues winding 
rightward through the channel indefinitely (Fig. 14c). 

Figure 14d illustrates all SZs (gray), SSU BFPs (red), 
and BIMs (red with arrows indicating the burning direc¬ 
tion) in this system. Two of the SZs contain the elliptic 
advective fixed points in the vortex centers. The others 
contain hyperbolic advective fixed points on the channel 
walls. Note that the SZs are slightly offset from a square 
lattice. This is due to the small wind added. We show 
only the SSU BFPs since, as we will show (Prop. 1), they 
are the generators of the unstable manifolds which com¬ 
bine to form FFs. Each BFP lies on a SZ boundary and, 
because it is SSU, is oriented away from the SZ. The 
BIMs emanating from these BFPs are similarly oriented. 
Each BIM spirals into a vortex center where it forms a 
cusp on an elliptic SZ (i.e. a SZ that contains an elliptic 
advective fixed point). Only the incoming portion of the 
cusp is pictured because, as shown earlier, the FF cannot 
contain cusps, and so the remainder of the BIM will not 
be relevant. 

Figure 14e summarizes the dynamical structures rele¬ 
vant to the behavior observed in Fig. 14b and 14c. The 
BIM core shown is responsible for blocking front propaga¬ 
tion at the bottom and center of the channel in Fig. 14b. 
The transverse stability of the BIM leads to the front’s 
rapid convergence upon it (Fig. 14b). As the front evolves 
further (Fig. 14c), it reaches the cusp at the end of the 
BIM core and winds around it. The BIM core does not 
form a complete span across the channel, and thus does 
not form a global barrier to the propagation of fronts. 
This is the situation seen in experimental images Figs. 4a, 
b, c. 

Now we increase the wind speed until it precisely bal¬ 
ances the burning speed, v w = vq (Fig. 15). Stimulat- 



FIG. 14. (Color online.) Small wind speed (r»o — 0.3, = 

0.15). (a) Fluid flow streamlines, fixed points and attached 
invariant manifolds, (b) Sequence of fronts shows preliminary 
convergence near bottom to roughly vertical curve, (c) Fur¬ 
ther evolution; lower edge converges to curved line while the 
rest proceeds around and to the right, (d) BIMs attached to 
BFPs. SZs shaded gray, (e) The one BIM most important for 
above front evolution - shown against advective structure. 


ing in the lower left (purple) we find that the reaction 
approaches a vertical curve (Fig. 15b), and so the reac¬ 
tion is confined to the left side. This appears to be a 
candidate for a frozen domain. In Fig. 15c, we test the 
stability of this region by introducing a small sinusoidal 
perturbation. The rightward component of this pertur¬ 
bation grows, eventually filling in the entire cell to the 
right, demonstrating that this region is not stable and 
therefore not a frozen domain. 

Let us examine the dynamical structures in Fig. 15d,e. 
The increase in wind has caused the SZs to shift slightly 
relative to Fig. 14d,e; the two on the lower boundary 
move together, as do the two on the upper boundary; 
those in the vortices move up or down depending on their 
rotational sense. The central BIM is now a straight ver¬ 
tical line. It is important to note that this BIM spans the 
entire channel with no cusp thus creating a global bar¬ 
rier to front propagation. Symmetry of the flow indicates 
that this BIM terminates at an SUU BFP at the top of 
the channel. We have previously argued that such a fixed 
point could not be on a FF, and it is this fixed point that 
leads to the instability demonstrated in Fig. 15c. 

Now we increase the wind beyond the critical value. In 
Fig. 16a a stimulation on the left converges to a burned 
region bounded by a smooth curve spanning the chan¬ 
nel. Unlike in Fig. 15c, a small rightward perturbation in 
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FIG. 15. (Color online.) Critical wind speed (vo = v w = 0.3). 
(a) Advective structure; similar to previous case, (b) This 
time, front progress (from the left) is completely blocked, (c) 
Perturbation of burned region shows instability, (d) Several 
BIMs, BFPs and SZs. (e) BIM of interest is a straight vertical 
line - coincides with separatrix of non-windy flow. 


Fig. 16b converges back to this smooth curve, and hence 
this curve is a FF. Figure 16c shows that the smooth 
bounding curve is the BIM emanating from the bottom 
BFP. Note that this BIM terminates at a point on the 
boundary that is not a BFP. This explains the situation 
seen in experimental images Fig. 4d, e, f as well as Fig. 5a. 

Now that we have seen BIMs act as both local and 
global barriers, we would like to understand the tran¬ 
sition between these two cases in more detail. Imag¬ 
ine a deformation that takes the BIM in Fig. 14e to 
the BIM in Fig. 16c; What might this deformation look 
like? Lemma 8 ensures that the angle between BIMs and 
streamlines is nonzero throughout the interior of the FZ. 
Therefore a BIM cannot form a tangency with the chan¬ 
nel wall (which must coincide with a streamline) in the 
interior of the FZ. Note, however, that a BIM cusp , on 
the boundary of a SZ, may encounter the channel wall 
without forming a tangency. In fact, this occurs when 
the cusp is perpendicular to the channel wall (Lemma 7). 
This observation suggests two deformation strategies: ei¬ 
ther move the existing cusp on the elliptic SZ to the wall, 
or create a new cusp on the hyperbolic SZ and slide the 
cusp to the wall. While the first mechanism seems more 
straightforward, and has not been ruled out theoretically, 
it has not yet been observed. However, the second mech¬ 
anism is observed here. 

In Fig. 17 we increase v w through the critical value 
v w = vq and follow the transformation of the BIM. Be- 



FIG. 16. (Color online.) Wind greater than critical (vq = 
0.3, v-u; = 0.4). A stimulation on the left (a) converges onto a 
smooth curve that spans the channel. In (b) a sinusoidal per¬ 
turbation of this curve converges back to the curve, implying 
that it is stable. (Only the last front is filled), (c) The BIM 
responsible for the FF spans the channel with no cusps. 


ginning with a subcritical v w value in Fig. 17a, we see 
the BIM (green) that comes up from BFP A (not shown) 
on the bottom wall and veers off to the right to form a 
cusp on the elliptic SZ (not shown). This cusp marks the 
end of the BIM core. 

Increasing the wind, the BIM is “blown backward” de¬ 
veloping a tangency (red and blue dashed) with the upper 
SZ. This tangency is not forbidden, because the SZ is not 
defined by a streamline. Since the front is burning away 
from the SZ, the tangency must occur at either an SSU 
or SSS BFP on the upper SZ (according to Lemma 6 and 
Fig. 25.) Because the SZ is convex in this case, the BFP 
must be SSU. The heteroclinic connection is illustrated 
by the coincidence of the unstable BIM of BFP A and 
the stable BIM of BFP B (red and blue dashed). 

Continuing to increase the wind, the BIM is blown 
further backward. Now it does not meet the SZ at a tan¬ 
gency, and so the heteroclinic connection is broken, giv¬ 
ing way to a cusp, the other option allowed by Lemma 6. 
This cusp slides along the SZ, with its angle changing to 
remain perpendicular to the fluid flow (Lemma 6). The 
cusp must rotate counterclockwise, at least initially, so 
that its tangent points into the SZ, as rotating clockwise 
would require the BIM to enter the SZ. 

The BIM soon arrives at another tangency with the 
SZ (Fig. 17b). Here, however, the BIM is burning into, 
rather than away from, the SZ. This tangency implies a 
heteroclinic connection with the SUU BFP C. (Again, 
the UUU case cannot occur because the SZ is convex; 
see Fig. 25.) In a generic three-dimensional dynamical 
system, heteroclinic connections between SSU and SUU 
fixed points are codimension-two occurrences. In this 
sytem, however, the BIMs are constrained to the two- 
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FIG. 17. (Color online.) The basic mechanism in the tran¬ 
sition to the first FF. We increase the wind speed, showing 
the interplay between the BIM from BFP A (not shown) and 
the upper SZ and its BFPs. (a) An SSU BFP B lies on the 
bottom of the SZ. Attached to it are BIMs (red) going left 
and right, both of which end in cusps on elliptic SZs (not 
shown), (i) A BIM (green) comes up from the SSU BFP 
A below (not shown) and then shadows the unstable BIM 
(red) going to the right, (ii) The BIM (red and blue dashed) 
forms a tangency/heteroclinic connection with the BFP B. 
(iii) The BIM (purple) is blown behind the heteroclinic con¬ 
nection, forming a cusp, (b) An SUU BFP C is shown at 
the top of the channel. A BIM (red) lies within its unstable 
manifold and goes off to the left, (iv) A BIM (green) slides 
leftward along the SZ, approaching the BFP C. (v) The BIM 
(red and blue dashed) forms a second tangency/heteroclinic 
connection with BFP C. (vi) The BIM (purple) is blown be¬ 
yond this heteroclinic connection, forming a complete span 
across the channel. Since the BFP positions and SZs change 
slightly with v w , the specific BFPs B and C shown, as well 
as their SZs, are calculated for the parameter values of the 
heteroclinic connections. 


dimensional sliding surface, and so the heteroclinic con¬ 
nection is a codimension-one occurrence. Said loosely, if 
a BIM is to sweep from one side of a SZ to the other, the 
BIM cannot avoid connecting with at least two BFPs on 
the SZ boundary. 

Increasing the wind still further, the BIM, blown en¬ 
tirely clear of the SZ, spans the entire channel, uninter¬ 
rupted by cusps. We have now arrived at the FF configu¬ 
ration in Fig. 16c. This FF topology persists for a signif¬ 
icant range of wind values. As seen in Fig. 18, the shape 
of this front can be nearly straight, or more boomerang¬ 
shaped, depending on the applied wind. Note that it is 
only due to the symmetry of the flow that the second 
heteroclinic connection in Fig. 17b occurs exactly when 
the BIM core first spans the channel 

At approximately wind value v w = 0.34 = 1.7i?o, the 
BIM encounters the upper left SZ in Fig. 18. Just like 
the BIM/SZ transition in Fig. 17, we must form a tan- 



FIG. 18. (Color online.) A series of FFs for increasing wind 
values (i>o = 0.3, Vo < v w < 1.7vo, blue to red). The BIM 
attached to BFP A is swept backward until it intersects BFP 
D. The SZs also shift and are colored accordingly. 


gency/heteroclinic connection (Lemma 6 and Fig. 25) to 
a new SSU BFP denoted BFP D. (Note the latter is 
rotated by roughly 7t/2 CCW in comparison to Fig. 17.) 
Foretelling this tangency, the red curves in Fig. 18 be¬ 
gin to curve upward near the upper channel wall. Once 
again, symmetry of the flow requires that BFP D be on 
the upper channel wall. 

As vq is increased still further, the BIM forms a 
cusp just behind the unstable BIM attached to BFP D 
(Fig. 19b), as seen in the mechanism in Fig. 17a. Note 
that while a front may wrap around the newly formed 
cusp attempting to bypass the initial BIM, it will shortly 
encounter the BIM attached to BFP D which has closed 
off this pathway (Fig. 19a). Here we have a FF that is 
composed of two distinct BIMs. Note that the burning 
region has a concave corner, with opening angle given 
by Lemma 8. The appearance of this concave corner 
is exactly what was observed in the experimental FF 
(Fig. 5b). 

The sequence in Fig. 20a takes the BIM through a se¬ 
ries of encounters with SZs as the wind speed is increased. 
(Here we consider Vo = 0.2 for simplicity.) The first en¬ 
counter is the attachment/detachment mechanism with 
the upper right SZ, analogous to that in Fig. 17. Here, 
however, the BIM detachment does not result in a BIM 
core that spans the channel. Rather, the BIM continues 
for some distance and then spirals in toward the ellip¬ 
tic SZ in the upper left, where it forms a cusp. As the 
wind speed increases, the cusp slides clockwise around 
the SZ until the BIM forms a new tangency with the up¬ 
per left SZ. The cusp on the elliptic SZ is “cut off” by 
this tangency, which dynamically precedes it along the 
BIM. This begins the mechanism of Fig. 17 again, after 
which the BIM core forms a complete span and defines a 
FF. 

Figure 20b shows a similar sequence as Fig. 20a for 
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FIG. 19. (Color online.) A composite FF formed from two 
BIMs. (i>o = 0.3,1^ = 0.525). The evolving front rapidly 
converges to BIM A and then winds around its cusp. How¬ 
ever, it is prevented from going further rightward by the short 
segment of BIM D. 


vq = 0.3. The main difference between these two images 
is that three SZs have merged into one in Fig. 20b. Con¬ 
sequently, the initial detachment of the BIM from the 
upper right of the SZ results in a new cusp formed near 
the bottom of the same SZ. Furthermore, as the cusp 
moves clockwise around the SZ, it is never “cut off”, but 
instead slides along the SZ to the channel wall. 

By flip-shift symmetry of the flow, the BIM attached 
to BFP D has undergone the same transition as the BIM 
attached to BFP A and so forms a FF as well (Fig. 21a). 
Importantly these two FFs intersect. Consequently, in 
addition to the frozen domains defined by single BIMs, 
the union of two neighboring frozen domains defines a 
distinct frozen domain. This union is continuously re¬ 
lated to the frozen domain observed in Fig. 19a as v w is 
increased. In Figs. 16 and 19, there is a 1-to-l correspon¬ 
dence between frozen domains and vortices in the chan¬ 
nel. Now in Fig. 21, the diversity of frozen domains (at 
fixed vo and v w ) has increased. We can have either a FF 
formed by a single BIM core (Fig. 21b), or by two inter¬ 
secting BIM cores (Fig. 21c, d,e). Note that the diversity 
of FFs in Figs. 21(b-e) is produced by small changes to 
the initial stimulation point. 

As the wind is increased, the process in Fig. 20 is re¬ 
peated. The BIM slides along the upper channel wall un¬ 
til it encounters an SSU BFP on a SZ. It moves around 
the SZ until it moves completely to the left of the SZ 
and reconnects to the channel wall. This process occurs 
once for each vortex pair. After each such occurrence, 



FIG. 20. (Color online.) Increasing the wind beyond first in¬ 
stance of FF generates new transitions, (a) (r»o = 0.2, 0.3 < 
v w < 0.6) Blue FF rapidly attaches to and detaches from 
SZ. Upon detaching, the BIM “jumps” straight to a cusp 
on an elliptic SZ. With even higher wind, another attach¬ 
ment / detachment leads again to a complete span, (b) 
(r»o s= 0.3,0.45 < v w < 0.8) Illustration of similar transitions 
where “jumping” occurs all within a single connected SZ. 

the BIM acquires a new intersection with another BIM 
emanating from the opposite side of the channel. We 
can thereby enumerate all frozen domains of this system 
for given values of Vo and v w . (Fig. 22). Finally, when 
min(|u|) > vq, there are no SZs, and therefore no BFPs, 
and therefore no FFs. 
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FIG. 21. (Color online.) FF diversity and sensitivity to ini¬ 
tial stimulation. (ao = 0.3, = 0.95). (a) BIMs A and D 

(related by flip-shift symmetry) each form a complete span, 
and intersect. Nearby stimulations (small pink dots near the 
right side) lead to different asymptotic frozen domains. The 
frozen domains fall into two classes: (b) and (c-e). (b) The 
FF is composed of a single BIM, which spans the channel, 
(c-e) The FF is composed of two BIMs. 



FIG. 22. (Color online.) For the windy alternating vortex 
chain flow we can enumerate the increasing number of possible 
frozen domains that occur with increasing wind speed. In this 
example, there are four FF shapes (up to flip-shift symmetry). 
(v 0 0.2, v w 1.18) 

VI. CONCLUSIONS 

The ability of a heterogeneous flow to freeze reaction 
fronts in the presence of an imposed wind appears to 
be quite general. Frozen fronts (“sustained patterns”) 
have been observed numerically in simulations of oceanic 
plankton blooms [35]. Frozen fronts have also been seen 
both experimentally and numerically in reacting flows in 
a porous media with a through-flow [36, 37] . We have also 


conducted experiments on frozen fronts in extended flows 
composed of two-dimensional arrays of vortices, either 
ordered or disordered [38]. As is the case in this paper, 
the frozen fronts in an extended flow with a wind are due 
to patterns of overlapping BIMs. 

This work suggests several directions of future re¬ 
search. In the context of design and control, this analysis 
could be used to develop a reacting fluid flow with some 
desired property. An obvious example is a system with 
maximal reaction rate. Given some class of accessible 
fluid flows, the reaction rate can be readily maximized 
by computing the lengths of FFs. Another example is 
reaction rate stability. We might be given a particular 
flow perturbation and seek the base flow that minimizes 
reaction rate fluctuation. 

It may be desirable to generate a FF with a partic¬ 
ular geometry. For instance, there may be a region in 
the neighborhood of the FF that we wish to keep strictly 
separated from the front (e.g. a sensor in the vicinity of 
a combustion front that cannot withstand the tempera¬ 
tures of the front itself). The analysis here provides a 
detailed connection between the stream function and FF 
shape making these questions accessible. 
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Appendix A: Sliding fronts 

Although we study fronts propagating in time- 
independent fluid flows, the fronts themselves certainly 
need not be time-invariant. For instance, a fast- 
propagating front in a weak flow will evolve approxi¬ 
mately as a circle of increasing radius. Loosely speaking, 
this is because each front element in the circle “burns 
beyond itself”. For a front to be time-invariant, each el¬ 
ement must instead “slide along itself”. In this section, 
we make this statement clear and derive several conse¬ 
quences. 

Definition 1. A front element, i.e. a point in xyO-space, 
is said to be sliding when r oc g, where g = [cos#, sin#]. 
Equivalently, 

r n = 0, (Al) 

where n = [sin #, — cos #]. 

While the sliding property is defined for any fluid flow, 
it is of most use when the flow is time-independent, as we 
have assumed throughout this paper and its appendices. 

The “sliding” constraint Eq. (Al) is illustrated geo¬ 
metrically in Fig. 9. For a given xy location, there are 



































13 




FIG. 23. (Color online.) Hyperbolic fluid flow, x = —Tx, y = 
-\-Ay. (v 0 = 0.35, A = 1) (a) Sliding surface, (b) Streamlines 
of the w + field, (c) Streamlines of the advective fluid flow, 
(d) Streamlines of the w_ field. 


FIG. 24. (Color online.) Elliptic fluid flow, x = — Ay, y = 
+Ax. (vo = 0.35, A = 1) (a) Sliding surface, (b) Streamlines 
of the w + field, (c) Streamlines of the advective fluid flow, 
(d) Streamlines of the w_ field. 


either zero, one or two solutions for 6 satisfying this con¬ 
straint. Where the fluid speed is small (|u| < i? 0 ), there 
is no solution; we call such a region a “slow zone” (SZ). 

Lemma 1. There are no sliding elements inside a SZ. 

Proof. Combining the sliding constraint Eq. (Al) with 
Eq. (2), we find |u • n| = |^o|. This cannot be satisfied 
for |u| < Vq- □ 

Where the fluid speed is large (|u| > vq), there are 
two solutions to Eq. (Al). Where the fluid speed equals 
the burning speed (|u| = vq), these two solutions are 
degenerate. We call a region where |u| > vo a “fast 
zone” (FZ). The sliding constraint Eq. (Al) defines a two- 
dimensional submanifold of xyO- space, called the sliding 
surface , which can be viewed as a double-branched sur¬ 
face over the FZ. Figures 23a and 24a show the slid¬ 
ing surface for a hyperbolic and an elliptic flow, respec¬ 
tively [39]. In Figs. 23a and 24a we see that, when viewed 
from above, these sliding surfaces have a hole in the mid¬ 
dle exactly where the SZ is. 

Lemma 2. If at any time some portion of a SZ is burned, 
the asymptotic burned domain will include that entire SZ. 


Lemma 3. Sliding is an invariant property. That is, if a 
front element is sliding, every element along its trajectory 
under Eq. (2) is also sliding. Hence we may speak of 
sliding trajectories. 

Proof. We examine the time derivative of Eq. (Al). 



= (uijrj + v 0 0gi)hi A- r k g k 6 
= UiUijfj + r k g k 0 


= ±|r |hiUijgj ± \r\g k g k 0 


= ±|r|(-0)±|r|0 = O, 


where the first equality follows from Eq. (2a) and the fact 
that dn/dt = g 0 , the second from the orthogonality of n 
and g, the third from the sliding assumption r = d=|r|g, 
and the fourth from Eq. (2b). □ 


A consequence of this lemma is that the sliding surface 
is dynamically invariant. 

Recall that a front is a curve (r(A), 0(A)) that every¬ 
where satisfies the front compatibility criterion Eq. (1), 
which is expressed equivalently as 


Proof. Within the SZ, the velocity of the front is every¬ 
where greater than the fluid. Therefore, no direction of 
motion is forbidden to the front, and so the front will 
eventually access all parts of the SZ. □ 


dr/dX • n = 0. (A2) 

Lemma 4. A trajectory of Eq. (2) is sliding if and only 
if the curve it sweeps out is a front. 
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Proof. Choosing A = t, Eq. (Al) is equivalent to 
Eq. (A2). □ 

In light of Lemma 4, we may refer to a sliding tra¬ 
jectory as a sliding front More generally, we make the 
following definition. 

Definition 2. A sliding front is a smooth curve that 
everywhere satisfies Eq. (Al), or equivalently Eq. (A2). 

Note that a sliding front may be composed of multiple 
trajectories joined at fixed points. Also, any segment of 
a sliding front is also referred to as a sliding front. 

Lemma 5. BIMs are sliding fronts, and thus lie within 
the sliding surface. 

Proof. A BIM is the unstable invariant manifold of an 
SSU BFP. Since we consider time-independent flows, this 
invariant manifold is also a trajectory. As shown in 
Ref. [8], BIMs satisfy the front compatibility criterion. 
Therefore, by Lemma 4, BIMs are sliding fronts. □ 

Since the sliding surface is invariant, it is natural to 
restrict Eq. (2) to this surface. We next derive an ex¬ 
plicit expression for this 2D flow. Applying the slid¬ 
ing constraint Eq. (Al) to the front element dynamics 
Eq. (2), we have u • n = — v$. Using this to resolve the 
unit vector n into components, we have n • u = — ^o/|u| 
and n • u x = =F\A — C L ’o/l u l) 2 , where rU or u 1 - is a 
righthanded rotation by7r/2 of uor u = u/|u|. Inserting 
the resolved form of n into Eq. (2) we have, 



(A3) 

This defines two flows over the FZ, one for each of the 
two branches of the sliding surface. Our sign convention 
is such that w + is the flow on the branch where r = +|r|g 
(the + branch), and w_ is the flow on the branch where 
r = —|r|g (the — branch). 

Equation (A3) shows that the w± fields are unde¬ 
fined (complex-valued) within the SZ and are zero on 
its boundary, confirming Lemma 1. All fixed points of 
Eq. (2) are fixed points of Eq. (A3) because all BFPs triv¬ 
ially satisfy the sliding constraint. These BFPs he on the 
SZ boundary. Equation (A3) also has a set of spurious 
fixed points at all other points along the SZ boundary, 
i.e. where |u| = vq. However, we ignore these spuri¬ 
ous fixed points as they are not physically relevant fixed 
points of Eq. (2), but rather result from the square-root 
singularity in Eq. (A3), obtained by projecting Eq. (2) 
onto xy- space. This square-root singularity also invali¬ 
dates the uniqueness of solutions to Eq. (A3) at the SZ 
boundary. Thus, there are other physically relevant tra¬ 
jectories that pass through the SZ boundary. 

Figures 23b, d and 24b, d illustrate the w± flows for the 
cases of hyperbolic and elliptic fluid flow, respectively. 
The gray regions are the SZs. BFPs are indicated on the 


boundary of the hyperbolic SZ. The stable and unsta¬ 
ble manifolds of these BFPs are shown in blue and red, 
respectively. 

Reference [8] proved that for any BFP, g is an eigen¬ 
vector of Uij, i.e. 


II 

is 

(A4) 

where g is the eigenvalue. Reference [8] also defined the 
quantity g !, 

g' = g + vqK, 

(A5) 

where ft is the signed curvature of the SZ boundary at 
the BFP. (ft < 0 means that n points toward the cen¬ 
ter of curvature.) Below, we reproduce Theorem 4 from 
Ref. [8], [40] 

Theorem 1. For a time-independent, 
flow u, the eigenvalues about a BFP are 

incompressible 

Ao = — Ab 

(A6) 

A± = - (-/X ± vV + 4/x/y 

) , (A7) 


where g and g! are given by Eqs. (A4) and (A5). The 
linear stability of a BFP is thus determined by the signs 
of g and g! according to the following table. 



g > 0 g < 0 

in! > 0 

sus uuu 

fj,' < 0 

sss suu 


We next specialize this result to the dynamics on the 
sliding surface. 

Corollary 1. The eigenvalues for a BFP of the dynam¬ 
ics Eq. (2) restricted to the sliding surface are given by 
X± from Eq. (Al). The xy-projection of each of the cor¬ 
responding eigenvectors is proportional to g. 

The sliding surface stability information is summarized 
in Fig. 25. For each of the four stability types, the first 
two stabilities (in bold) describe the dynamics within the 
invariant sliding surface. Equation A5 places restrictions 
on the local convexity of the SZ at the BFP. These pos¬ 
sibilities are illustrated in Fig. 25. 

Lemma 6. At an intersection p between a sliding front 
and the boundary of a SZ, the fluid flow is perpendicular 
to the sliding front. If the sliding front is tangent to the 
boundary, p is a BFP. Otherwise, p is a cusp along the 
sliding trajectory. All BFPs and cusps of sliding fronts 
occur at the intersection between a sliding front and SZ 
boundary. 

Proof. Combining the sliding constraint Eq. (Al) with 
the front element dynamics Eq. (2) gives u n = —vq. 
Since |u| = Vo on the SZ boundary, u = — t> 0 h. Thus u 
is perpendicular to g. 

In Thm. 2 of Ref. [8], it was shown that a necessary 
and sufficient condition for a BFP was for it to be on the 
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FIG. 25. (Color online.) Four BFP stability types. Black ar¬ 
row indicates burning direction. Gray regions are SZs. Two 
un/stable manifolds of each BFP are within the sliding sur¬ 
face. Each manifold is labeled with +/— indicating its corre¬ 
sponding branch of the sliding surface. 


boundary of the SZ with n perpendicular to the bound¬ 
ary. Thus g tangent to the boundary implies a BFP. 

If g is not tangent to the boundary, then since the slid¬ 
ing trajectory cannot enter the SZ, it reaches the bound¬ 
ary and then must reverse direction forming a cusp. 

Finally, a BFP and cusp both require r = 0. This 
satisfies the sliding constraint and also implies |u| = v$. 

□ 

Lemma 7. A BFP or cusp on the boundary of the fluid 
domain, i.e. at a wall, must have g perpendicular to that 
boundary. 

Proof. This follows from the previous Lemma 6 and that 
the fluid velocity of an incompressible fluid is tangent to 
the fluid domain boundary. □ 

An incompressible 2D fluid flow can be specified by a 
stream function 4/(r), with u x = d^/dy,u y = —d^/dx. 
Each fluid element follows a level set, or streamline, of T. 
Front elements, on the other hand, do not follow stream¬ 
lines, but generally cross them one way or the other de¬ 
pending on their relative orientation. 

Lemma 8. Sliding fronts cross streamlines such that 
#n,u, the angle between n and u ; satisfies cos(#n jU ) = 

-V|u|. 



FIG. 26. (Color online.) Sliding fronts represented on the 
graph of the stream function T near: (a) a hyperbolic point, 
(b) an elliptic point. Sliding fronts are curves of constant 
increase (or decrease) in the value of the stream function. 


crossed. To this end, we calculate the rate at which 4/ 
changes when viewed from the frame of an individual 


front element. 


DT _ d T dx d T dy d T 

Dt dx dt + dy dt + dt 

(A8) 

= - u y x + u x y 

(A9) 

= —(y + v o cos 6)x + (x — Vo sin 0)y 

(A10) 

= —u 0 r • g 

(All) 

= — sgn(r • g)u 0 r , 

(A12) 


where the last equality makes use of the sliding con¬ 
straint. Scaling by the front element speed, we find the 
simple relation 


— = - sgn(r • g)v 0 , (A13) 

where s measures the euclidean x?/-length along the tra¬ 
jectory. Equation (A13) shows that the sliding trajecto¬ 
ries on the + (—) branch of the sliding surface are those 
that climb the stream function with the constant rate of 
descent (ascent) vq. It is straightforward to show that the 
only front elements that ascend or descend at a constant 
rate are sliding. 

Consider two streamlines with values 4/ a and \lT and 
a sliding front F a $ that connects one to the other with 
no intervening cusps. The xy-length of the segment F a ^ 
follows directly from Eq. (A13), 


Proof. From Eq. (2), u = r — non. The sliding front 
condition implies u = =L|r|g—non. Dotting with n, n-u = 
—n 0 . □ 

Lemma 8 means that sliding fronts are never tangent 
to streamlines, except in the no/|u| -T 0 limit. Physical 
boundaries of the fluid (channel walls) are particularly 
important streamlines at which this lemma can be uti¬ 
lized. 

In addition to the angle at which sliding fronts cross 
streamlines, we can examine how rapidly they are 


\Fa,b\ 



(A14) 


This expression is particularly useful when thinking 
about FFs in channel flows of arbitrary geometry. Since 
the channel wall enforces a boundary condition of con¬ 
stant 4/, the FF length is found through Eq. (A14). Thus 
the length of a FF that spans a channel depends only 
on this “energy difference” between the two walls and 
the burning speed, and not on other details of the flow. 
While Eq. (A14) was derived for a single sliding front, 
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it also holds for FFs that are composed of multiple BIM 
cores. This can be seen by applying Eq. (A14) to each 
BIM segment separately. Interestingly, this implies that 
multiple FFs existing in the same flow must have the 
same length, even in the absence of any flow symmetry. 

Equation (A14) also implies that a channel of width 
W cannot support a FF if AT < vqW. Furthermore, if 
we assume that a flow u without wind gives no net flow 
down the channel, then for the flow u + v w , Eq. (A14) 
becomes 


\F a ,b\ = Wv w /v 0 . (A15) 

This can be interpreted as the equality of fluid flux across 
the FF and across the channel width. Eq. (A15) shows 
that v w > vo is necessary but not sufficient. In the case 
that FFs do occur at v w = t?o, they must be straight 
lines that meet the channel walls at right angles. This 
occurs exactly when the original fluid flow (v w = 0) has 
a vertical advective separatrix. This condition is met 
for the windy alternating vortex chain model; addition¬ 
ally Ref. [10] experimentally demonstrated that v w = vo 
marked the onset of FFs. However, it is not difficult to 
construct flows where v w = vo is not sufficient for the 
existence of FFs. 


Appendix B: Stability of Frozen Domains 

At the beginning of Sec. IV, we specified that frozen do¬ 
mains should be stable under small perturbations. Here, 
we define this stability more precisely. For a given in¬ 
variant burned domain, with boundary F, we define an 
allowable distortion of F at a point r E F to be a dis¬ 
tortion such that F remains unchanged outside a ball of 
radius e(r) > 0 centered at r. Note that the value of e(r) 
is not fixed but may vary with the point r. 

Definition 3 (Stability of frozen domains/fronts). A 
frozen domain/front is required to be stable in the follow¬ 
ing sense. There must exist a function e(r) of each point 
r along the front F, i.e. the domain boundary, such that 
after any allowable distortion, the front remains point- 
wise close to F and converges pointwise to F as time 
goes to infinity. 

Here pointwise close is in the “Lyapunov” sense, in 
that the maximum (over all time) distance from the time- 
evolved distorted front to the original front F remains 
bounded and goes to zero as e goes to zero. 

With this definition, one can easily verify the argu¬ 
ment of Sec. IV proving that a FF cannot contain an 
SUU or UUU BFP. One also sees that a FF can contain 
an SSU or SSS BFP, and that the curves constructed in 
Prop. 1 are stable. Regarding the latter, it is interesting 
to note how a perturbation of the FF F returns to F. 
First, a perturbation localized to the neighborhood of an 
SSS point simply shrinks in size, back into the original 
FF, due to the SSS point’s being a sink. Consider now 



FIG. 27. (Color online.) The fates of localized perturbations 
to a FF. a) The perturbation eventually shrinks into an SSS 
BFP on the FF. b) The perturbation strikes a wall of the fluid 
domain, c) The perturbation strikes a second BIM core at a 
concave corner of the FF. d) The perturbation goes to infinity 
along a BIM core that stretches to infinity. Though the size of 
the perturbation could grow indefinitely, the perturbed front 
still returns to the FF pointwise. 


a perturbation localized at some point r of the FF that 
is not an SSS BFP. This perturbation will be “swept” 
along the front, away from the unstable BFP that gener¬ 
ates the BIM; the perturbation might even initially grow 
in size. The localized perturbation will continue to follow 
the BIM segment from which it was perturbed, and will 
subsequently encounter either an SSS point, a domain 
wall, another BIM segment of the FF, or it will be swept 
to infinity (Fig. 27). In the initial three cases, it is clear 
that the perturbation will disappear as it either shrinks 
into the SSS BFP (Fig. 27a), strikes the wall (Fig. 27b), 
or runs into the already burned region (Fig. 27c), assum¬ 
ing the initial size of the perturbation e(r) is sufficiently 
small. The case in which the perturbation is swept to 
infinity (Fig. 27d) requires an additional assumption on 
the far field nature of the fluid flow, addressed in Sec. C 4. 


Appendix C: Dynamics at infinity 

In this paper, we assume that the fluid velocity field 
u(r) is smooth and either (i) has a bounded stream func¬ 
tion or (ii) is “localized” with a simple far-held behavior. 

More precisely, case (i) assumes that u(r) is generated 
by a stream function T(r) that has a global maximum 
and minimum over the fluid domain. This case applies, 
for example, to the stream function in Eq. (4), since in 
our restricted domain y is bounded between 0 and 1. Un¬ 
der case (i), no sliding front may be infinitely long with¬ 
out striking a SZ, since a sliding front element increases 
or decreases its stream function value at the constant rate 
uo, Eq. (A13). Thus, for case (i) there are no BIM cores 
that stretch all the way to infinity. 

The remainder of this appendix concerns case (ii), 
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stated more precisely as the requirement that, as r goes 
to infinity, u(r) behaves as a homogeneous polynomial 
P k (x,y) of power k > 0 in the variables (x,y), i.e. 

u(r)=P k (x,y) + E(x,y), (Cl) 

where E (x,y) is a term that grows more slowly than r k . 
More precisely, we require that 

hm r _/c E(r cos </>, r sin</>) = 0, (C2) 

r—> oo 

with the polar angle <p held constant in the limit. Fur¬ 
thermore, we require that radial derivatives of E grow 
more slowly, according to 

r\£ 

lim r _/c+ ^ 7 —rE(r cos0, r sin0) = 0, 0 < L (C3) 

r—>• oo QT^ 

This assumption on the far field behavior eliminates, for 
example, infinite arrays of vortices [though such arrays 
could be allowed under case (i)], while allowing many 
important cases, including flows with constant far field 
velocities, linear hyperbolic and elliptic flows, other poly¬ 
nomial flows, and flows constructed from arbitrary con¬ 
figurations of a finite number of vortices. 


1. Existence criteria for fixed points at infinity 

By a burning fixed point “at infinity”, we intuitively 
mean a front element trajectory that attains a constant 
orientation 6 and polar angle <p at an infinite value of r. 
We can formalize this definition by adapting the standard 
Poincare compactification of the plane [41]. First, map 
the radial distance r, 0 < r < oo, to a new radial variable 
p, 0 < p < 7r/2, defined by 

p = arctanr. (C4) 

Similarly, introduce scaled Cartesian coordinates p = 
(pxipy) — ( x ,y)p/ r - This transformation maps the xy- 
plane to an open disk of radius 7r/2 in the p-plane. We 
then include in our phase space the boundary of the disk 
at p = 7r/2, representing the “circle at infinity”. Equa¬ 
tion (C4) yields 

p =- J — (I - p® p)r + (cos 2 p){p® p)r (C5) 

tanp 

= y — (I -p<8>p)(u + v 0 n) (C6) 

tan p 

+ (cos 2 p)(p 0 p)(u + u 0 n), 

where the second equality follows from Eq. (2a). Here, 
I is the identity matrix and p 0 p is the tensor product 
with components (p 0 p)ij = PiPj/ p 2 . Equation (C6) is 
smooth in (p x ,p y ,0) everywhere except possibly at the 
boundary p = n/2. 

Note that Eq. (Cl) can be reexpressed in (p, p) coor¬ 
dinates as 


~ k (p)Q(cf>) + E (p,4>), (C7) 


where Q is smooth in p and where 


lim 

p—>-7r/2 


cos^ £ (p) 


d^_ 

dp £ 


E(p, p) = 0, 


o < e, 


(C8) 


with fixed p. Thus as p goes to 7r/2, Eq. (C6) scales 
as 1/cos /c_1 p, and is thereby singular for k > 1. To 
remove this singularity, the ODE time-parameter t can 
be replaced by a new scaled parameter s defined by 


ds 

dt 


cos^ -1 p’ 


k > 1. 


When k = 0, we make no scaling. This yields 

d pcos k p, t ^ 

— p =—-(I - p0p)(u + u o n) 

ds sin p 


(C9) 


(CIO) 


+ (eos /c+1 p)(p 0 p)(u + u 0 n), k > 1. 


The right-hand side of Eq. (CIO) is now smooth over 
the entire closed disk p < tt/2 and all 0 values, except 
possibly at p = 0, which does not concern us since our 
focus is at p = 7r/2. One can also see that 

4# = _( cos fc-i p)fi i u i jg j , k> 1, (Cll) 


is smooth in (p x ,p y ,0) for all p < 7r/2, but might be 
singular at p = tt/2. 

We define a fixed point at infinity of Eq. (2) to be a 
fixed point (p*,p*,6>*) at p* = tt/2 for Eqs. (CIO) and 
(Cll) (k > 1) or Eqs. (C5) and (2b) (k = 0). We next 
restrict attention to those fixed points at infinity that 
are the limits of sliding fronts extending to infinity. The 
tangent direction of a sliding front that converges upon a 
fixed point at infinity must thus converge to d= p. Thus, 
such a fixed point must yield 


g* = ±p, 

•2 


lim —- 

r—^oo r r 


= 1 . 


(C12) 

(C13) 


Furthermore, notice that at p = 7r/2, the term 

(cos^ -1 p)uijCjj in Eq. (Cll) will be noninfinite if g = 4=p 
[see Eqs. (C7) and (C8)], a fortuitous consequence of 
searching for fixed points that are the limits of sliding 
fronts. 

Equation (C13) can be rewritten as 


. ( u r + v 0 n r ) 2 

1 = hm —-—— = hm 


oo |u + Uon|- 


= lim 


r ^°° |u zb Vo4>\ 2 
cos - 2k (p)Q 2 r {4i) 


(C14) 


p-Hr/2 COS ~ 2k (p)Q%((f>) + [COS ~ k (p)Qct,((f>) ± V 0 ] 2 

(C15) 

Q 2 M) ( C16 ) 


= lim 


P-S-tt/2 Q 2 (<^) + [Q </,(<!>) ± VO COS fc (/3)] 2 ’ 


where the second equality follows from Eq. (C12) and the 
third from Eqs. (C7) and (C8). This implies that for a 


u(p, <p) = cos 
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fixed point at infinity (that is the limit of a sliding front) 
in polar coordinates (tt/2, 0*, 6*) 

Q*{<n = ±«o cos fc (7T/2) = |^° * * = J’ (C17) 

Assume now that k > 1. We seek necessary and suf¬ 
ficient conditions for dp/ds = 0 at p = tt/2. Equa¬ 
tions (C7), (C8), and (CIO) imply that at p = tt/2 

^p(tt/2, 4>, 0) = |(l - p 0 p)Q(^) = 

(C18) 

which is independent of 0. Thus, a necessary and suffi¬ 
cient condition for dp/ds = 0 at p = 7r/2 is 


*) = 0, fc>l, 


(C19) 


consistent with the sliding front condition Eq. (Cl 7). 

We now seek necessary and sufficient conditions for 
d6/ds = 0. Based on Eq. (Cll), we define the scaled 
Jacobian matrix 


J(p, 0 )= COS /C p 




(C20) 


expressed in the (r, 0) basis. Under assumptions (C7) 
and (C8), one can show that J takes on the following 
form at p = tt/2. 


J(tt/2,0) 


2. Stability of fixed points at infinity 

Assume k > 1. To determine the stability of fixed 
points at infinity, we compute the Jacobian matrix of the 
flow (CIO) and (Cll) in (p, 0, 6) coordinates, i.e. 


J = 


d(p, <t>, 0) 


(C25) 


Applying the results of Theorem 2, it can be shown that 
the Jacobian at a fixed point at infinity (to which a sliding 
front converges) is 


J = 


Qr 


0 


0 kQ^^ Q(f) : cf) kQ r 


where 


0'p = —k(cos p) k 1 (E ( p + Von#). 


(C26) 


(C27) 


Applying the area-preservation constraint V • u = 0 to 
Eqs. (C7) and (C8), we find 


— (1 T k^Q r , 


from which 





1 

1 

O 

0 

kQr^ffl ^5r,(/)(0) 

(C21) 

J = 

0' p —(1 + k)Q r 

0 

kQ 0 ( 0 ) ^5</>,</>(0) 


0 — (1 T k^jkQ r 

— (1 + 2k)Q r 


(C28) 


(C29) 


For a fixed point at p = 7r/2, we apply Eq. (C19) to Thus, the eigenvalues of J all have the same sign, which 
obtain is the opposite sign to Q r . We thus have the following. 


J(tt/2,0*) 


kQ r {r) QrA </>*) 

o Q^m 


(C22) 


Setting d6/ds = 0 in Eq. (Cll) is equivalent to g being 
an eigenvector of J. From Eqs. (Cl2) and (C22), how¬ 
ever, we see that g = p is already guaranteed to be an 
eigenvector of J(7r/2,0*). Thus, d9/ds = 0 is already 
implied by Eqs. (C12) and (C19). 

Assume now that k = 0. Then P k (x,y) is constant, 
and it is easy to see from Eq. (C6) and Eq. (2b) that 
p = 0 and 6 = 0 at p = 7r/2 for any 0 and 6. However, 0 
and 0 are no longer arbitrary when we apply the sliding 
front conditions Eq. (C12) and Eq. (C17). 

To summarize, we have the following conditions on the 
existence of fixed points at infinity. 

Theorem 2. Assume the incompressible flow u satisfies 
Eqs. (Cl) and (C3) or equivalently Eqs. (C7) and (C8). 
Then a sliding front ends in a fixed point at infinity with 
orientation and position anqles 0 * and 0* if and only if 
g(6>*) = ±p(<p*) and 


QM) = 0, for k> 1, 
Q<p{4>*) = ±uo> for k = 0. 


(C23) 

(C24) 


Theorem 3. Assume the incompressible flow u satis¬ 
fies Eqs. (Cl) and (C3), with k > 1. Suppose a sliding 
front ends in a fixed point {tt/2 , 0 *,#*) at infinity and 
assume Q r (0*) ^ 0. Then the fixed point has stability 
SSS (UUU) if the flow along the sliding front is radially 
outward (inward), i.e. Q r (0*) > 0 (Q r (<!>*) < 0). 

In the k = 0 case, the eigenvalues of the Jacobian 
are all 0, and stability must be determined by nonlinear 
analysis. 


3. Behavior of frozen fronts at infinity 

Consider a sliding front that extends to infinity with¬ 
out containing a cusp beyond some set radius. Such a 
front must be a trajectory of the w± field Eq. (A3) with 
a single choice of sign. Transforming this vector field into 
the p-coordinates according to Eq. (C5) and scaling time 
via Eq. (C9), we obtain a smooth 2D vector field over 
p -space in the neighborhood of the boundary p = n/2. 
Thus the sliding front that extends to infinity, i.e. con¬ 
verges to p = 7r/2 in the p-coordinates, can do only one 
of two things. If the boundary p = 7r/2 contains no 













19 


fixed point, then the boundary is a limit cycle to which 
the sliding front converges. Otherwise, the sliding front 
must converge upon a fixed point at p = tt/ 2. 

Let us suppose that the boundary p = 7r/2 is a limit 
cycle, which we also suppose is stable. (The analysis of 
the unstable case is similar.) Then, a circle atp = 7r/2 —e 
will converge outward to 7r/2 for a small enough e. In the 
original ^-coordinates, this means that a large enough 
circle of front elements would expand outward without 
bound. If these front elements were passive tracers of the 
flow, this outward expansion would clearly violate conser¬ 
vation of area of the underlying flow. Similarly, if these 
front elements were all pointing inward with nonzero Vq, 
they would clearly still violate conservation of area of 
the underlying flow. Thus, the front elements can only 
be facing outward if the circle is ultimately to grow in 
size. In summary, we have the following. 

Proposition 2. Assume the incompressible flow u sat¬ 
isfies Eqs. (Cl) and (CS). A sliding front that moves 
outward to (inward from) infinity, without cusps beyond 
a certain radius, with its burning direction pointed inward 
(outward) must converge upon a fixed point (7r/2,0*,#*) 
at infinity. 

Consider now a frozen front that converges on a limit 
cycle at infinity. There would in fact need to be two 
separate frozen fronts converging upon a limit cycle at 
infinity, one burning inward and one burning outward, 
to create a burned strip that spirals outward to infinity. 
However, the previous proposition shows that only one 
such burning direction is possible. Hence, we have shown 
the following. 

Proposition 3. Assume the incompressible flow u sat¬ 
isfies Eqs. (Cl) and (CS). Any frozen front that extends 
to infinity, without any concave corners or BFPs beyond 
a certain radius, must be a sliding front converging upon 
a fixed point at infinity. 


4. Stability of frozen fronts extending to infinity 

We now complete the discussion of frozen front stabil¬ 
ity begun in Appendix B. In particular we consider a FF 
F that extends to infinity and an allowable perturbation 
to F at a point r G F that is, at least initially, swept 
toward infinity. Our concern is that this perturbation 
may not only grow in time, but could burn far enough 
away from F that it would no longer be carried away to 
infinity but would leave some new part of the fluid, not 
in the original burned domain defined by F, burned for¬ 
ever. This problem is resolved, however, by the Poincare 
compactification. If the sliding front converges on a sta¬ 
ble fixed point at p — 7r/2, then one can always find a 
sufficiently small perturbation e(r) such that the pertur¬ 
bation never grows too large to be collapsed by the SSS 
BFP at infinity. 



FIG. 28. (Color online.) When the parameter c is below 2, 
the curvature of the SZ is not strong enough to create the 
SSU to SSU-SSS-SSU bifurcation. ( Vo = 1,/x = 1 ,c = 7/4). 
(a) SSU BFP with attached sliding fronts. Three other BFPs 
on the SZ boundary are not shown, (b) Fronts on the sliding 
surface seen in 3D from the center of the SZ. 


Appendix D: SSS points 

It was shown in Ref. [8] that SSS is a possible stability 
type for BFPs. Here we provide an explicit example for 
the fluid flow surrounding an SSS point. 

The stream function 


-r c/i 2 s 

* = -voy - wy + -z—y 

6v 0 


(Dl) 


produces the following flow field with a BFP at (x,y) = 
(0,0) and g = y 


2 

c/r 2 

U x = -v 0 -fj,x+ —y , 
2v 0 

u y = (iy, 


(D2) 

(D3) 


where vq is the burning speed, p is the eigenvalue from 
Eq. (A4), and c is a new parameter. Note that this 
flow can be understood as the combination of a uniform 
“wind”, a linear hyperbolic flow, and a Poiseuille flow. 
The stability of the BFP can be determined from Theo¬ 
rem 1 once p! is computed. Ref. [8], Eqs. (19) and (31), 
showed that 


m' = (2m 2 -I’oa-gV/x, (D4) 

a- g = nigjgkUijk, (D5) 

from which it is straightforward to compute a g = cp 2 /v o 
and p! = (2 — c)p. Thus, the BFP at the origin has 
stability SSS if p > 0 and c > 2. If in addition c < 9/4, 
Eqs. (A6) and (A7) show that the three eigenvalues are 
real. 

Figure 28 illustrates this flow for c = 7/4, which is 
too small to produce an SSS BFP at the origin. Equa¬ 
tion (A5) allows us to interpret this fact as due to the 
magnitude of the SZ curvature \k\ being too small. Topo¬ 
logically, the dynamics has the same structure as the lin¬ 
ear hyperbolic flow Fig. 23, since the origin is an SSU 
BFP. The BIM (red) forms a FF extending to infinity. 

Increasing c to 8.5/4, the magnitude of the SZ curva¬ 
ture at the BFP is increased (Fig. 29). The SSU BFP in 
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FIG. 29. (Color online.) (i>o = 1,/jl = l,c = 8.5/4). (a) 
When the central SSS BFP has real eigenvalues, the BIMs 
emanating from the two neighboring SSU BFPs reach it with 
no cusps, (b) BIMs run very close to the sliding surface fold. 


FIG. 30. (Color online.) (vo = 1,/jl — 1 , c = 15/4). (a) When 
the central SSS BFP has complex eigenvalues, the BIMs em¬ 
anating from the SSU BFPs repeatedly overshoot the SSS 
BFP, and form cusps, (b) From side view, it is clear that 
these cusps are caused by the local spiraling of trajectories. 


Fig. 28 has bifurcated into an SSS BFP at the origin sur¬ 


rounded by two SSU BFPs. The BIMs emanating from 
the SSU BFPs toward the central SSS BFP actually ter¬ 
minate on the SSS BFP (most easily seen in Fig. 29b.) 
Thus, all three BFPs lie on the FF composed of these 
BIMs. 

Increasing c further to 15/4, the SSS point acquires 
complex eigenvalues. Therefore, any sliding trajectory 
that reaches the center point must first encircle it an 
infinite number of times (Fig. 30). Because the SSS point 
is on the fold of the sliding surface (where |u| = uo), 
such a trajectory must pass through the fold an infinite 
number of times, generating an infinite sequence of cusps. 
Therefore, for any point along a sliding trajectory (which 
we imagine to be a part of some FF), there must exist an 
infinite number of cusps between it and the SSS point. 
This ensures that an SSS point with complex eigenvalues 
cannot lie on a FF. 

Conversely, an SSS point can only lie on a FF if its 
eigenvalues are real. Such SSS points do exist in the flow 
Eq. (D2) for parameters satisfying vo > 0,/i > 0,8/4 < 
c < 9/4. Thus in Fig. 29 the real eigenvalues ensure that 
the SSS point lies on the FF attached to the neighbor¬ 
ing SSU points. Because the eigenvalues are real for an 
open interval in parameter space, the connection between 
the SSU and SSS BFPs in Fig. 29 is structurally stable. 
Here the FF runs very close to the SZ boundary, however 
there is a small gap. Constrast this with the SSU to SSU 
connection previously shown in Fig. 13, which was not 
structurally stable. 
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